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Abstract 



Dimerization of a spin-half Heisenberg antiferromagnet on a square lattice is investigated by taking unexpanded 
exchange couplings. Several dimerized configurations are considered some of which are shown to have lower ground 
state energies than others. In particular, the lattice deformations resulting in alternate strong and weak couplings 
along both the principal axes of a square lattice are shown to result in a larger gain in energy. In addition, a 'columnar' 
configuration is shown to have a lower ground state energy and a faster increase in the energy gap parameter than 
a 'staggered' configuration. The inclusion of unexpanded exchange coupling leads to a power law behaviour for the 
magnetic energy gain and the spin-Peierls gap, which is qualitatively different from that reported earlier. Instead 
of varying as S x , the two quantities are shown proportional to 8 V / |ln <5| . It is thus proposed that the logarithmic 
correction, that was regarded as an outcome of including umklapp processes for small distortions, can also be a direct 
consequence of using an un-truncated spin-spin exchange interaction. The calculations, which employ the coupled 
cluster method, lead to a conclusion result that the dimerization of a spin-half Heisenberg antiferromagnet on a sqaure 
lattice is unconditional. 



PACS numbers: 75.10.Jm, 74.65. +n, 75.50.Ee 
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I. INTRODUCTION 



It is known that dimerization lowers the ground state energy of a spin-half isotropic Heisenberg antiferromagnet 
|p]-^0|. In other words, the system stands to gain energy by such lattice deformations that render it dimerized with 
alternate weaker and stronger bonds between up and down spins on neighboring sites. On the other hand the lattice 
distortions cost energy and it is the net energy balance that would determine whether the gain in magnetic energy is 
large enough to affect the spin-Peierls transition through dimerization. In a phenomenological theory, this is usually 
seen in terms of an exponent showing the dependence of magnetic and elastic energies on a parameter 5, describing 
the extent of lattice deformations. It is related to the displacement of the ith atom through m = \(— 1) L S. The 
spin-dimer formation is usually described by the Hamiltonian 

H = j£[l + (-l) i (J]S i -S i+1 (1) 

i 

with alternate stronger and weaker bonds J(l + S) and J(l — S). 

Since < S < 1, and since elastic energies go typically as 5 2 , therefore if the magnetic energy gain varies with 5 
with an exponent less than 2 then in the limit 5 — > 0, the gain would overwhelm the cost and a spontaneous and 
unconditional dimerization should occur. If, however, the exponent is equal to or greater than 2 then dimerization 
may occur, conditional to the details of the interaction parameters. 

Such aspects as these have been studied extensively in Heisenberg antifcrromagnetic chains in which the critical 
exponent has been shown to favor an unconditional spin-Peierls transition, as summarized in Table 1. This aspect 
has also been revealed by experiments on quasi-one dimensional Heisenberg antiferromagnet CuGeO 



Table 1 lists the scaling laws proposed for the magnetic energy gain in one- dimensional antiferromagnets both 
in the small 5 (near critical) regime as well as the far critical regime (0 < S < 1). Of particular interest to us is 
the logarithmic correction for small (5, which is reported to be a result of including umklapp processes [pi. We find, 
however, that if, instead of the approximated form of the spin-spin exchange coupling, J(l ± <5), its full form 

^«)=|4 (2) 

is used then the magnetic energy gain comes to vary with <5 as tj^-jp exactly as calculated by the various methods 
shown in Table 1. Additionally, it is easy to see that this scaling behaviour will not be confined to the near critical 
region, but will operate over the entire range of S. The full form of the interaction includes the often-used truncated 
forms in an obvious way: when the distance between a pair of spins decreases from a to a(l — d), the exchange coupling 
is taken to increase from J to approximately J(l + 5), etc. In what follows, we shall use for exchange interaction the 
form in Eq.(^|). 

The situation in two-dimensions is a little more involved because of the additional factor of frustration due to a 
competing antiferromagnetic second neighbour interaction which can in principle destroy any LRO of the Neel type as 
well as the possibility of dimerization. Much of the study of two-dimensional Heisenberg antiferromagnet has therefore 
remained focused on the destruction of order by frustration. 

The matter of frustration aside, a simple dimerization of a square lattice is interesting in its own right because the 
lattice distortions can take place in more than one way, each one of the possible configurations giving a different de- 
pendence of the ground state energy on the dimerization parameter. Studies of dimerization of a spin-half Heisenberg 
antiferromagnet on a square lattice have mostly considered a columnar configuration, illustrated in Fig. 1(a) [^4|-^9|. 
That is to say, the lattice deformations and the consequent spin-singlet pair formations are taken to occur only along 
one of the square axes, say the x-axis. The nearest neighbour distance along that axis is taken to vary alternately as 
a(l — S) and a(l + 5), while that between neighbours in the perpendicular (y-) direction remains a. It has been shown 
that, like the chains, this configuration also gives rise to a continuous reduction of ground state energy with S. 

One can think of other possible configurations of two-dimensional dimerization which are different from this one in 
the sense that they allow for lattice distortions or changed spin-spin couplings along the y-direction also. We shall 
propose a few such configurations in the following, three of which are shown in Figs.l(b)-(d). 
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In Fig.(b), the lattice deformations along the x-direction is alternated as in Fig. (a), but the sequence of alternations 
is itself alternated as one goes along the y-direction. To distinguish it from the columnar dimerization, we call it stag- 
gered dimerization p3] . It makes the exchange coupling along the y-direction also dependent upon the dimerization 
parameter S. While the coupling along the x-direction is alternately and y^y, it is uniformly — ^=== along the 

y-direction. 

In contrast to the configurations (a) and (b), those in Figs. 1(c) and (d) allow for simultaneous dimerization along 
both x- and y-directions in the plane. The difference between (c) and (d) is the same as that between (a) and (b): 
configuration (c) is columnar and (d) is staggered. The former is called plaquette configuration |T5| , pO| |. These four 
configurations of a dimerized square lattice consisting of N spins are therefore characterized by the following nearest 
neighbour interactions. 



Configuration (a) 

J ^ = TTTxsJ^ J(1-XS), A = ±l 

Jy (/. 



That is to say, the dimerized configuration is described by the Hamiltonian 

VN r 



h = jj: 

Configuration (b) 

J y = vi+s 2 — J ( l _ t) 

and the Hamiltonian is given by 



h =jy: 

i,3 



1 



Si, 7 ■ Sj+i j + Si j ■ Si 



3 + 1 



(i + (-iy+is) 



Sij ■ Si + ij 



1 



Sij • Sijj-j-l 



Configuration (c) 



T — •' 
Jx > x ~ (1 + A<5) 



J(l - XS), A = ±1 



■ / «.a= (T^y- J(i + xs) 

with the Hamiltonian 



h =jy: 



(i + (-i)^) SiJ ' Sl+hl + (i + (-iys) SlJ ' S * J+1 



(3) 



(4) 



(5) 



Configuration (d) 



^,A = TT ^y-^(l-A ( 5), A = ±l 



and the Hamiltonian 



J (l - XS - (1 - 4r)<5 2 



' jv 



1,3 



(l + i-iy+'S) h3 



St+i 



5 2 + (1 + (— 1)J"<5) 2 



Si.j ' Sjj_|_i 



(6) 



3 



We would like to investigate the four configurations in order to see (i) which one of these leads to the largest gain 
in magnetic energy as the dimerization sets in, (ii) whether such a dimerization is conditional or otherwise, and (iii) 
if the use of un-truncated exchange coupling leads to a single scaling law valid for the entire range of S. 

Being a prototype of the superconducting phase transition plf , the appropriate order parameter for spin-Peierls 
transition must also be the energy gap between the ground and excited states - the energy required to break a singlet 
pair. We shall therefore look at the critical behaviour of the system under dimerization in terms of the (5-dependence 
of the energy gap parameter. 

A number of methods can be chosen for this purpose. Spin wave theory, either modified through Takahashi con- 
straint of zero magnetization or a Hartree-Fock approximated non-linear theory, is known to give surprisingly good 
results for spin-half Heisenberg antiferromagnet. Or, a spin wave theory in the spinless fermionic representation 
through Jordan- Wigner transformations takes care of fermionic correlations among the s — | spins. Coupled cluster 
method has also been extensively, and successfully, used for spin-half Heisenberg antiferromagnet in one and two space 
dimensions. 

The first two methods belong to the class of mean field theories and hence are not expected to be very reliable when 
it comes to determining critical exponents. The coupled cluster method, on the other hand, is a perturbation method 
in which increasingly higher order correlations can, in principle, be incorporated at will, and which has been shown 
to give satisfactory results even in the lower orders of perturbation. We believe that the coupled cluster method must 
be sufficiently good to see if the energetics allow a spin-Peierls transition to the alternative configurations proposed 
here. 



II. APPLICATION OF THE COUPLED-CLUSTER METHOD 



In the coupled cluster method it is first necessary to define a ket state starting from a model state | (f> >, which in 
our case is the Neel state. The exact ground state | \& > of the system can then be postulated as 

| * >= e s | > (7) 
where S is the correlation operator defined for an N particle system as 

S = E 5 « ( g ) 

71 

withS„ = Si 1 ,....,i»« Cj n (9) 

ii •■ -in 

and C\ is the creation operator defined with respect to the model state. The ground state energy can then be found 
as the eigenvalue of the Hamiltonian in the proposed ground state 

He s | >= E g e s \ 4> >. 

Taking inner product with < (j) \ e~ s gives 

E g =<(f>\ e~ s He s \ (f> > . 

The product e~ He can be written as a series of nested commutators in the well-known expansion 

e- s He s =H+[H,S]+j ] [[H,S],S} + (10) 

where in the present case the series terminates after the fourth term. 

It is usually easier to deal with the s = \ Heisenberg Hamiltonian by applying a rotation of 180° to the up spin 
sublattice such that in it S x — > — S x , S y — > S y and S z — » — S z . The fermionic character of the spin operators is then 
preserved by expressing them in terms of Pauli matrices: = \a 3 , j = x,y,z [ p2| , ^3[| . A general expression for the 
nearest neighbour spin Hamiltonian in 2D is then 
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where p is a vector to the four nearest neighbours. Correspondingly, the string operator S n can now be defined as 

S ^ = (iE E s =: i - j,..o„4< • • • <^ • • • (i 2 ) 

ii....i„ ji....j„ 

where subscripts i and j distinguish between sites on the two sublattices. We note that for spin half (<7j + ) 2 = 
(af) 2 = 0. Truncat ion of the summation up to the desired level gives rise to different schemes of approximation. 
Taking interaction only between the spins on adjacent sites gives the so-called SUB2-2 scheme. Including interactions 
with the second and fourth neighboring sites gives what is termed as SUB 2 _4 scheme. And taking the previous two 
schemes including interaction among the four adjacent sites give us what has been termed as local SUB4, or LSUB4 
for short. Each one of these approximations accounts for a different order of perturbation calculation, and takes into 
account a different order of inter-particle correlations. It has been noted that LSUB4 is a sufficiently good approxi- 
mation for calculating the ground state properties of a spin-half Heisenberg system [p3[ . 

Consider a general case: a Hamiltonian which has four different coupling constants for nearest neighbour interactions 
in two space dimensions. It can be written as 

^N/2 

H = —~ E/ ^ X ' X a 2i,j ' a 2i+\,j + Jy,X &i,2j ■ <?i,2,j+\] ■ (13) 

i,j A=±l 

Here i and j are the two components of the site indices on a square lattice. The correlation operators in the LSUB4 
scheme are defined as 



S 2 
S3 



i 



/ n a 2i+v,j + 9 n a ti- v ,i + h i[ a t]+v + 1 n 



v=0 



(14) 



In these equations, the coefficients ai^i^etc, are various forms of the coefficient Si x ; in the expressions for 

S2n- The ground state energy within the LSUB4 approximation comes out to be 



E„ 







Jx, + 1 


(i +o 





-'-Jx.-i {- + b 1 ) + J V!+1 (-+ci) + J a _i i^+dx) 



(15) 



The coefficients ai, a 2 , ■ • •, I are obtained as solutions of a set of coupled nonlinear equations. These equations arise 



from the fact that such matrix elements as < 



O 



>He £ 



> are all zero when the operator O is any product 



of creation operators, particularly if it is one of the operator products in the correlation operator S above. 

'He s >= 0; 
• 3 He s >= 0; 



< a i,2j a i,2j+u e SHpS >= °; 



< (T 2ij (T 2i+v,^ 



v = ±1,±3 
v = ±1,±3 
v = ±1 
v = ±\ 



'i,2j u i,2j+v u i,2j+2v u i.]+3v^ ") v — (16) 

where S = S 2 + S3 + S4 . These equations translate into the following twelve equations for the unknown parameters: 
(a{ - 1 - 2a 3 b 1 - 2f)J x>+1 + (2oi + 2oi&i - 2a 3 6 3 )^,-i = 
(26! + 20161 - 2a 3 6 3 ) J Xj+1 + (6 2 - 1 - 2 fll 6 3 - 2g)J x ,- t = 
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(2a 3 + 2aia 3 ) J x . +1 + (2a 3 - a\ + 2a 3 bi - f)J x -i= 
(26 3 - b\ + 2 ai b 3 - g)J x , +1 + (2b 3 + 2b 1 b 3 )J x ,- 1 = 

(-2a 3 6i + 2ai/ + a 3 &f - ai<z 3 6i)J X)+ i+(/ - a\ + 2h\j + a 3 .g + 2aia 3 6 3 ) J x _i = 
(# - 6f + 2ai.g + 6 3 / + 2a 3 &i& 3 ) J^+i + (-2ai& 3 + 2b^ + a\b 3 - ai&i& 3 ) J x -1 = 
(c? - 1 - 2c 3 di - 2h)J y . +1 + (2ci + 2c x di - 2c 3 d 3 ) J y _i = 
(2d! + 2c 1 d 1 - 2c 3 d 3 )J y , +1 + (d? - 1 - 2cid 3 - 2Z) J y _i = 
(2c 3 + 2cic 3 ) J Vt+ i + (2c 3 - cf + 2c 3 di - h)J y -i = 
(2d 3 - d\ + 2 Cl d 3 - l)J y , +1 + {2d 3 + 2d 1 d 3 )J y ,- 1 = 

(— 2c 3 di + 2c\h + c 3 d\ — cie 3 di) J^+i + (h — cf + 2d x h + c 3 l + 2ac 3 d 3 )J y -i = 
(l-dj + 2c x l + d 3 h + 2c 3 d 1 d 3 )J y , +1 + (-2cid 3 + 2d x l + c\d 3 - cieM 3 ) J y _i = 

Setting all the coupling constants J M equal reduces the number of equations from twelve to three and yields exactly 
the same equations as obtained by others |2^,|2^]. The two sets of six equations each independently determines the 
six coefficients contained in each of them. As expected, the equations are symmetric in some coefficients. The twelve 
coefficients are to be evaluated by solving the above coupled equations numerically for each of the configurations 
separately by substituting appropriate values of J Xy \ and J Vi \. 

To calculate the energy gap we shall construct the excited ket state |\I/ e > in term of a linear excitation operator 
X, which, operating on the ground state >i takes the system to an excited state: >= X >= Xe \<j> > . 
This operator is constructed as a linear combination of products of creation operators p3| 

X = J>„ (17) 

n 

with 

x n = £ xj,...j„44-<- ( 18 ) 

The first excited state is obtained by the operator 

j 

where j can be any site of the two sublattices. It is easily seen that the first excitation energy is 

E e = i ( ~ + 01 + 61 + a s + b 3 ^j {J x .+i + J x _i) . (20) 

The energy gap for a given 5 is A(<5) = E e (S) — \E g (8)\ ■ The order parameter for the spin-Peierls transition is therefore 
D{5) = A(<5) — A(0). This is the energy required to break a dimerized singlet pair. 



I. III. RESULTS 



The ground state energy can now be calculated as a function of the dimerization parameter 8 . Previous calculations 
have invariably taken spin-spin exchange couplings alternately as J(l±5), which, as mentioned above, is an expansion 
of the interaction in Eq.(|J) to order 6, implying that the results are valid only in the limit <5 — > 0. We notice in our 
calculations that if all the expansions are terminated beyond the order 8 then the distinction between configurations 
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(a) and (b) [Eqs.(|3|, ^)] disappears. On the other hand, if the expansion is taken to one order higher, then there 
remains no way to distinguish between configurations (c) and (d), Eqs.(||,^). We must therefore either go to orders 
beyond S 2 in the expansion, or retain the interactions in their unexpanded form Eq.(||). We do the latter. An added 
advantage is that the results will then be valid in the limit 5 — > 1. 

Our calculations confirm that, like the chain, the ground state energy of all the four configurations of a dimerized 
spin-half Heisenberg antiferromagnet on a square lattice decreases with 5. This is shown in Figures 2, where e(d) — e(0) 
is plotted against 6 for the proposed configurations. The conclusion is not new for the simple columnar dimerization 
of Fig. 1(a) and (c) |b|-|l9|]. However, what is significant is that the ground state energy goes down with 6 more rapidly 
for configurations (c) and (d), Figs. l(c)-(d). In fact, Fig. 2 shows that the (5-dependence is markedly different for the 
two types of dimerized configurations: one in which dimerization takes place only along one axis, and the other, in 
which it occurs along both the directions. The rate of decrease is significantly higher for the latter. Also, for the entire 
range of S the columnar configurations lead to a greater gain in magnetic energy than the staggered dimerization. The 
two observations put together show that the plaquettc configuration of Fig. 1(c) is energetically the most favourable 
state, as noted earlier Particularly in the complete range of 5 (0 < S < 1), the plaquette configuration stands 
out as the most preferred one, while there is hardly a discernible difference among the other three. 

It is worth pointing out here that the much simpler mean field methods of spin wave theory - either in the bosonic 
representation through Holstcin-Primakoff transformations, or in the fermionic representation through Jordan- Wigner 
transformations - yield very similar results. This has been checked by us separately. 

To see if the spin-Peierls transition setting in the Heisenberg antiferromagnet on a square lattice is conditional or 
otherwise, we look at the balance between the gain in magnetic energy through dimerization and the cost in energy 
for elastic deformations; the latter going as S 2 . Earlier calculations on a square lattice had obtained the exponent as 1 
[|l8t for columnar configuration (a) , and as 2 for columnar as well as staggered configurations Jl5|,[]l| . Our results must 
be different from these because instead of J(l ±5), we take the unapproximated exchange coupling J(ry) = t^tT' 
However, they must agree with the earlier results in the limit of small 5. 

We find that the gain in magnetic energy does not quite scale with S as a simple power law: it scales as jj^-jyfor 
all the four configurations. The exponent v = 1.5 in the range < 8 < 0.1, and is equal to 1 in the complete range 
< S < 1. This is true for all the four configurations. That Tang and Hirsch had obtained a simple power law with 
exponent 2 is easy to understand if one notes that the difference between S 2 and tj^t is difficult to discern in the 
small 8 regime, particularly while fitting to a curve. 

Since j^-jy for v = 1 or 1.5 is always larger than d 2 in the entire range < S < 1, the spin-Peierls transition in a 
spin-half Heisenberg antiferromagnet on a square lattice must be unconditional. 

As shown in Table 1, the dimerization of an antiferromagnetic chain scales as yj^-jy, but only in the small S regime 
(the near critical regime). In chains, the factor of yj^y is taken to serve as a correction due to umklapp processes. 
This is not the case in a square lattice. In fact our CCM results show that this may not be the case even in chains 
when the exchange couplings in the dimerized state are taken as instead of the approximated J(l ± 5). We find 

for chains that the magnetic energy gain fits to t^ji m the entire range of <5 rather than only in the range of small 5. 

With the full exchange couplings, the coefficients come out for the chain to be v = 1 for < 5 < 1, and v = 1.3 —1.6 
for < S < 0.1. These give decent comparisons with the numbers in Table 1. 

The S dependence of the order parameter for the four configurations is shown in Fig. 3. An obvious conclusion from 
these results is that the dimerized state becomes more stabilized with increasing <5. Like the magnetic energy gain, 
the order parameter D also scales with S as yj^-jy with the same values for v in both the small as well as full region 
of S. 

The difference between the dimerization of a square lattice along only one direction (Fig. 1(a) and (b)) and along 
both the directions (Fig. 1(c) and (d)) is again markedly brought out in Fig. 3. Also the columnar configurations 
again appear to be the preferred modes of dimerization over the staggered configurations for having higher values of 
the order parameter in the region of small 6. 

To summarize, we have studied the spin-Peierls dimerization of a spin-half Heisenberg antiferromagnet on a square 
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lattice taking unapproximated exchange couplings starting from the ansatz J{rij) — We have included different 
possibilities of dimerization. The ground state energy decreases continuously with increasing dimerization for all the 
proposed configurations, . Of the four configurations, those with dimerization taking place simultaneously along both 
the principal square axes have markedly lower ground state energies than those with dimerization along only one 
of the axes. Also, those with columnar dimerization have consistently lower energies than those with the staggered 
dimerization. The spin-Peierls gap also corroborates the above conclusions. It has also been shown that the magnetic 
energy gain as well as the spin-Peierls gap scale with the dimerization parameter <5 as j^jj . It has also been asserted 
that this form of scaling does not need to be a consequence of corrections due to umklapp processes, and that the 
same can be said of the spin-Peierls transition in chains. As a result of this, it is concluded that the spin-Peierls 
transition in a spin-half antiferromagnet on a square lattice is unconditional. 

* E-mail: aiman@physics.sdnpk.undp.org 
nayyar ©physics . sdnpk . undp . org 
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FIGURE CAPTIONS 



Figure 1: Four configurations for the dimerization of a square lattice, (a) Columnar: the nearest neighbour coupling 
along the horizontal direction alternates between J(l — 5) and J(l + 5), while that along the vertical direction remains 
J. (b) staggered: like (a), the dimerization occurs along one direction only, but the sequence of alternate couplings 
itself alternates along the other direction, (c) Dimerization along both the directions, making a plaquette of four 
nearest neighbour spins, (d) Again dimerization along both the directions, but taken staggered along the vertical 
direction. 

Figure 2: The gain in magnetic energy e(S) — e(0) as dimerization sets in with increasing S for the four configura- 
tions; (a) in the range < 5 < 0.1, and (b) in < 5 < 1. 

Figure 3: Dependence of the energy gap parameter D on S for the four dimerization configurations; (a) in the range 
< S < 0.1, and (b) in < 6 < 1. 
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Table- 1 : Summary of the critical exponents for spin-Peierls transition in a Heisenberg chain determined by 
methods. 



Method 


Interval 


e(8) - e(0) 


Exponent 




< 8 < 1 




x — 4/3 


T? pnoTTTifi li 7i\ 1~i on oTmin 

JL LC1JU1 J.11CU JZjCI L1UJ1 (LlULlU 1— r 


< 8 < 1 




x = 1.53 


2-level RG § 


< 8 < 1 
0.05 < <S < 0.1 
0.4 < 8 < 0.5 


<5 2 7| ln(5) | 
<5 2 7| ln(<5) 


x = 1.78 
2i> = 1.68±g"|| 
2v = 1.31 ±0.02 


Excitation spectrum l@ 


< 8 < 1 


<P 


x = 1.36+^ 


Valence bond ID 


8 < 0.05 
8 > 0.05 


<5 2 7|ln( 

S x 


v = 2/3 
a; = 1.36±^i 


Finit size scaling 8 


< 8 < 0.1 
< <5 < 1 


8 2v /\ ln{8)\ 
S x 


v = 0.71 ±0.01 
x = 1.34 ±0.02 


Exact diagonahzation l§ 


< 8 < 0.1 


8 2v / | \n{8)\ 


v = 2/3 


DMRG 


.5 < 0.05 


8 X 


x = 2/3 
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